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Abstract 

I discuss the shape of the high energy end of the spectrum of particles arising from 
diffusive shock acceleration in the presence of (i) additional diffusive escape from the 
accelerator, (ii) continuous energy losses, (iii) energy changes arising from interactions. 
The form of the spectrum near cut-off is sensitive to these processes as well as to the 
momentum-dependence of the diffusion coefficients and the compression ratio, and so the 
spectrum of any radiation emitted by the accelerated particles may reflect the physical 
conditions of the acceleration region. Results presented in this paper have applications 
in interpreting the spectral energy distributions of many types of astrophysical object in- 
cluding supernova remnants (SNR), active galactic nuclei (AGN) and acceleration sources 
of ultra high energy cosmic rays (UHE CR). Except for extremely nearby sources, spectral 
features imprinted on the spectrum of UHE CR during the acceleration process will be 
largely eroded during propagation, but the spectrum of UHE neutrinos produced in inter- 
actions of UHE CR with radiation, both during cosmic ray acceleration and subsequent 
propagation through the cosmic microwave background radiation, contains sufficient in- 
formation to determine the cut-off momentum of the UHE CR just after acceleration for 
reasonable assumptions. Observation of these UHE neutrinos by the Pierre Auger Obser- 
vatory may help in identifying the sources of the highest energy cosmic rays. 

PACS Numbers: 98.70. Sa, 95.85.Ry, 98.38.Mz, 98.54.Gr 

Keywords: cosmic rays, acceleration, propagation, neutrinos, supernova remnants, radio 
galaxies 
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1 Introduction 



For particle acceleration by electric fields induced by the motion of magnetic fields B 
(including those at astrophysical shocks), the maximum rate of momentum gain by rela- 
tivistic particles of charge Ze can be written (in SI units) 



dp 
'di 



mZecB (1) 



where ^{p) < 1 and depends on the details of the acceleration mechanism. For example, 
in diffusive shock acceleration (DSA) ^{p) will depend on the shock velocity and diffusion 
coefficients. I shall refer to ^ (p) as the acceleration rate parameter. Detailed and rigorous 
treatments of DSA are given in several review articles |I1I21E1; see particularly the review 
by Jones & Ellison j3] , on the plasma physics of shock acceleration, which also includes a 
brief historical review and refers to early work. 

Protheroe and Stanev jH] proposed a box model of DSA which was able to reproduce 
the essential features of DSA in a simple leaky-box scenario. Particles were injected 
into the box at momentum pq, and while inside the box their energy increased at the 
average rate that would apply given the physical parameters of the shock acceleration 
being modelled, i.e. the upstream and downstream diffusion coefficients, the shock velocity 
and compression ratio. Particles being accelerated leak out of the box, again at the average 
rate that would apply given the physical parameters, and these escaping particles represent 
the accelerated particles escaping downstream in the standard shock acceleration picture. 
This simple scheme is very easy to implement in a numerical or Monte Carlo program for 
investigating the structure of pile-ups and cut-offs in the case of adding additional leakage 
terms (e.g. simulating the finite size of any accelerator), energy losses during acceleration 
(e.g. synchrotron loss) or particle interactions such as inverse Compton (IC) losses in the 
Klein-Nishina regime. These processes could be important in determining the detailed 
shape of spectrum of non-thermal radiation from SNR or AGN. 

Where protons are accelerated, interactions with magnetic fields, matter or radiation 
may determine where and how the spectrum cuts off. Of current interest is the hoped-for 
observation of vr*^ TeV gamma rays from SNR from pp interactions of cosmic rays freshly 
accelerated at SNR shocks. At higher energies, the UHE CR may interact with background 
radiation fields, including the cosmic microwave background radiation (CMBR), by pion 
photoproduction interactions both during acceleration and propagation. This leads to 
a predicted spectral cut-off, the Greisen-Zatsepin-Kuzmin (GZK) cut-off EI at about 
10^" eV for particles which have travelled more than a few tens of Mpc from their sources. 
However, several experiments have reported UHE CR events with energies above 10^° eV 
[S| with the highest energy event having 3 x 10^° eV jH), although the latest HiRes data 
[TU] seems to show a cut-off consistent with the expected GZK cut-off. 

Information about the acceleration environment is, in principle, imprinted in the shape 
of the high energy end of the spectrum of accelerated particles. In the case of the UHE CR, 
one of the main applications of the present work is to see to what extent this information 
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remains, either in the cosmic ray spectrum, or in the spectrum of secondary particles such 
as neutrinos produced during both acceleration and propagation [TTl |T2l UH] • Hence it 
is important to calculate these fluxes resulting from cascading in the CMBR. Spectral 
data at energies above lO^^eV and directional results, notably from the AGASA project, 
are very suggestive of fascinating, unexpected physics ^3]. A new era in the fleld is 
beginning with the commissioning of the Pierre Auger Observatory ^3] which comprises 
a pair of 3000 km^ arrays (one under construction in Mendoza Province, Argentina, and 
one planned for Utah, U.S.A.) employing both particle and optical detectors. The Auger 
Observatory will also have the capability of detecting UHE neutrinos which could carry 
signiflcant information about the sites of acceleration (see refs. \W\ IT7j for recent review 
of UHE CR). 

I shall flrst examine the case of continuous energy losses, then consider how the spectral 
features are altered for the case of interactions with various mean inelasticity values. As a 
speciflc example, I shall take the acceleration and propagation of the UHE CR, and briefly 
re-examine the likely acceleration sites of UHE CR. That will be followed by a description 
of Monte Carlo simulations of DSA using the box model and incorporating Bethe-Heitler 
pair production and pion photoproduction interactions in the CMBR; the results obtained 
being applicable also to other temperature blackbody flelds with suitable scaling. Finally, 

1 shall propagate these spectra to see whether any useful features remain in the UHE CR 
spectrum, and whether or not the spectrum of secondaries such as neutrinos can be used 
to infer physical conditions in the acceleration region of UHE CR. 

2 Roll-ofFs and cut-offs in the spectrum of acceler- 
ated particles 

The inclusion of escape processes, in addition to the usual escape downstream leads to 
a smooth roll-off in the spectrum of accelerated particles which is most pronounced at 
momenta near the momentum at which the two escape rates are equal, which I shall refer 
to as the "maximum momentum", Pmax- This smooth roll-off occurs over up to three 
decades in energy and the spectrum can extend beyond Pmax- Its shape depends on the 
momentum dependence of the diffusion coefficients, on the momentum dependence of the 
additional escape rate, and on the intrinsic differential spectral index, F, of the accelerated 
particles which would apply for the case of no energy losses or additional escape processes. 
Such roll-offs have very recently been noted in test-particle Monte Carlo simulations of 
shock acceleration at shocks in a cylindrical jet geometry where there is sideways leakage 
out of the jet [TH] . 

When continuous energy losses are included the spectrum is cut-off sharply at a mo- 
mentum at which the total rate of momentum change is zero, and I shall refer to this as 
the "cut-off momentum", pcut- Depending on the intrinsic spectral index, and momentum 
dependence of the diffusion coefficient, a pile-up may occur before the cut-off. Protheroe 
and Stanev j3] subtracted from the momentum-gain rate a term representing the energy- 
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loss rate. Drury et al. JH] showed that since the physical size of the "box" increases with 
energy, synchrotron losses can cause a particle in the downstream region to effectively 
"fall out" of the box, and that this process could be represented by an additional term in 
the formula for the escape rate. This term is included in the present work, and I confirm 
the result of Drury et al. jTH] that the pile-up is less, particularly for low compression 
ratios. 

In the box model of shock acceleration, particles of momentum po are injected into the 
acceleration zone, or "box", and while inside the box are accelerated at a rate r^cdp) and 
escape from the box at a rate resc{p)- These two rates uniquely determine the spectrum of 
accelerated particles, i.e. those escaping from the box. The rate of change of momentum 
during shock acceleration in the presence of energy losses is 



pri,cc{p) = P [rga.m{p) - rio 



dp 




dp 




dp 




'dt 
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(2) 



Note that, by definition, rgain(Pcut) = ?^ioss(Pcut)- 

The momentum gain and escape rates depends on the diffusion coefficients, which 
are usually assumed to have a power-law dependence on momentum, k{p) oc p^, with 
the exponent depending on the nature of the turbulence present in the magnetic field: 
5 = 1/3 (Kolmogorov spectrum), 1/2 (Kraichnan spectrum) or 1 (completely disordered 
field). When the diffusion coefficients upstream and downstream have the same power-law 
dependence on momentum, Ki(p) = Kip^ and K2{p) = K2P^, the momentum gain rate is 



gam 



(p) (X p , giving 



racc(p) = '^loss(Pc 



P 

yPcut, 



- noss(p)- 



(3) 



In the case of dp/dt oc losses, such as synchrotron loss, rioss(p) oc p and one has 



ip) '"loss {p cut I 



p 

yPcut, 



p 

Pcut 



(4) 



If there are no energy losses and no additional escape processes a power-law spectrum 
results, and the integral spectral index , (F-l), equals the ratio of the escape rate to 
the acceleration rate. In this case, the escape rate representing escape downstream is 
(r — l)r\oss{Pcnt)(yP/Pcnt)~^ ■ The acceleration zone extends distances Li{p) = Ki{p)/ui 
and 1^2 (p) = i^2{p)/u2 upstream and downstream from the shock, respectively. Adding 
an extra escape term r™^^(p,pmax) responsible for the "maximum momentum", and an 
escape term rlf^"^^{p) resulting from energy losses causing particles to fall out of the box 
gives the total escape rate 



'"esc(p) = (r - l)rioss(Pcut) 



P 

,Pcut, 



+ reT(P,Pmax)+r::r*(p) 



(5) 
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For the case of a constant extra escape rate as adopted by |S] 

Cr(P, Pmax) = (r - l)n,UPcnt) , (6) 

\ Pent J 

and for the case of diffusive escape out of the box, e.g. diffusion perpendicular to the 
shock normal to the edge of the acceleration region with diffusion coefficient k oc p^, 

CsTiP, Pmax) = (r - l)rioss(Pcut) ( ^] • (7) 

This latter case is probably more realistic than constant escape, and will give weaker 
roU-offs. 

The escape rate resulting from energy losses is 



jaiiout/..A 1 dp 



r'~iP) 



L{p) dp dt 



(8) 



where L{p) = Li{p) + L2{p), and only the decrease of L2{p) with decreasing momentum 
contributes to escape as particles falling out of the upstream part of the box (length Li{p)) 
will be advected back into the box. For Ki{p) = Kip^ and H2{p) = K2P^, 

r':^'"'\p)=Si2noss{p) (9) 
where £2 = L2{p)/L{p) is related to the compression ratio R = Ui/u2 = (2 + r)/(r — 1), 

K2R 2 + r 



K^ + K2R {KJK2){T ~l) + 2 + T' 



(10) 



Drury et al. jTH] appear to use K2 = Ki (possibly to represent a parallel shock for which 
in an ideal case B2=Bi), giving £2 = (2 + r)/(l + 2r) = 4/5 for a strong shock. For a very 
weak shock £2 — * 1/2. For a perpendicular shock, for which in an ideal case B2 = RBi, 
one might expect K2=Ki/ R, giving £2 = 1/2. Except where otherwise stated, I shall use 
£2 = 1/2 in this paper. 



3 Spectrum for cut-ofF due to continuous losses 



At time t after injecting Nq particles of momentum pq into the acceleration zone, the 
number of particles remaining in the acceleration zone, N{t), is obtained by solving 



dN 
'dt 



which has solution 

N{t) = No exp 



'"esc 

\p{t)]dt 



Nq exp 



- / resc{p)^dp 
Jpo dp 



'12) 



Hence, the spectrum of escaping particles is 



^ dp dt dp °pracc(p) 



ip) ' 
dp 



Po pr^ccip) 



(13) 



For continuous dp/dt oc losses (synchrotron loss or IC scattering in the Thomson 
regime), and defining x = p/pcnt, from Eq. UHlthe spectrum for the case of xq <^ x < 1 ^ 
Xmax is given by 

- ^^^^^S^ (t)"'"' - .-0'--*""- . (14, 

For continuous dp/dt oc p losses (e.g., adiabatic), from Eq. UHlthe spectrum for the case 
of Xmax > 1 is givcu by 

6i2X + (F - l)x-' fx\''- f x-' - 1 ^"-^)/^+^^ 



X^ ^ — X \XqJ \Xq — 1 

Note that in both cases, for x <C 1 this reduces to F^scip) ~ (r — ^)xo~^x~^ , as it must. 
For Xmax oo, and F = 1.5, 2.0 and 2.5 the shape of the cut-off is shown in Fig. [T] For 
synchrotron losses (Fig. [T^) the spectrum is hardly affected up to a; = 0.3 and strong 
pile-ups are present for F < 2, but for the case of adiabatic losses (Fig. [TJd) a pile-up 
occurs only for F < 2 and 6 = 1. Furthermore, adiabatic losses cause a much stronger 
turnover and affects the spectrum at significantly lower energies, down to x ~ 10~^ for 
6 = 1/3. Of course the situation is modified for finite Xmax- In that case, the integral in 
Eq. 1131 involves Hypergeometric Functions (as used in ref. j3]) and these are obtained by 
summation of an infinite series. 

The spectral index as a function of momentum is given by 

dlnFesc _ p dr esc _ P dr^^^c ^csc(p) ^-^^-^ 

dhip Tcscip) dp racc(p) dp r^cc{p)' 

One can obtain (the logarithm of) the differential energy spectrum by integrating Eq. IT^ 
numerically from Inpo to In p. For the case of inclusion of both synchrotron losses plus 
an extra momentum-dependent escape term (Eq. [7j), Eq. ITHl gives 

d In Fesc _ 6i2X + 6{T - l){x'x-li - x'^) 6x-' + x 
dlnx 6^2^ + (F — l)(x~'^ + x''x~ax) ~ ^ 

_ 1 _ 6i2X + iT-l){x-' + x'x-fJ ^ 

X~' — X 

and results for various Xmax, 6 and F values are shown in Fig. |21 These results apply 
equally to electrons and protons if synchrotron radiation is the dominant energy loss 
process. An application in which the cut-off in the proton energy spectrum is due to 
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synchrotron loss is given by the Synchrotron Proton Blazar model for the jets of BL Lac 
type active galactic nuclei which have weak accretion disks [201 1^- In contrast, in quasars 
with strong accretion disks, the spectrum of accelerated protons may be cut off by pion 
photoproduction interactions on thermal photons from the accretion disk. Cut-offs due 
to pion photoproduction will be discussed in Section 6. 

Finally, one can find the criterion for a pile-up by examining the spectral index as 
x — > 1 pni- The denominators of the 2nd and 4th terms of Eq. El tend to zero as a; — > 1, 
and so the spectral index becomes infinite. The sign of the denominator is always positive 
if S > —1, and so whether the spectral index tends to +oo or — oo depends on the sign of 
the numerator. Hence a strong pile-up occurs if 

5(1 - 4) + 2 - r - (r - i)x;Z > O- (is) 

Taking x^aax oo and ^2 = (2 + r)/(l + 2r), as assumed by Drury et al. [TH], I obtain 

(2-5) + (3 + 5)r-2r2 > (19) 

from which the maximum values of F for a pile-up are 2.069 {S = 1/3), 2.106 {S = 1/2), 
and 2.225 {6 = 1). Taking the compression ratio to be i? = (2 + F)/(F — 1), for 6 = 1 
a pile-up is expected if i? > 3.44949 (in agreement with Drury et al. (19j). For ^2=1/2, 
as adopted in the present paper, the conditions for pile-ups for synchrotron losses and 
adiabatic losses are F < 2 + 6/2 and F < 1 + 6, respectively, and are seen to be consistent 
with the spectra in Fig. ^ 

4 Spectrum for cut-ofF due to interactions 

Interactions of relativistic particles with matter or radiation generally result in the inter- 
acting particle losing energy, which is either given to the struck particle or photon (elastic 
collision) or used in production of secondary particles (inelastic collision). If the mean 
interaction length (mean free path) is Xmtip), for ultra-relativistic particles the rate of 
interaction is rint(p) = c/Aint(p)- If the mean inelasticity, i.e. average fraction of energy 
lost per interaction, is a{p) then the effective loss rate for these non-continuous losses is 
noss(p) = a(p)'^int(p)- Of course, if a{p) <C 1 then one can approximate the interactions as 
continuous energy losses. I shall now investigate how the spectrum of particles subject to 
interactions during acceleration differs from that of particles subject to continuous energy 
losses with the same effective loss rate. 

In the simulation with Nq particles injected, a particle is injected at time t=0 with 
momentum p(0) = Po = XoPcut and statistical weight w{0) = Wq = 1/Nq, and its sub- 
sequent momentum, p(t) = x{t)pcut, and weight, w(t), are determined after successive 
time steps At chosen to be much smaller than the smallest time-scale in the problem; 
At = 0.01xo'^/rioss(Pcut) is used in these simulations. In each time step, first the momen- 
tum is changed, p{t + At) = p(t)[l + Atri,cc{p)], and then the probability of escaping in 
time At is estimated as Pesc = {1 — exp[— Atresc(p)]} where resc here does not include the 
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'"esc""* term. Then w(t)Pesc particles with momenta p are binned in a histogram of accel- 
erated particles, and the particle's weight is changed to reflect the fraction not escaping, 
w{t + At) = w{t){l — Peso)- Next, the probability of interacting in time At is estimated, 
-Pint = {1 ~ exp[— At rint(p)]}, and a random number is generated to determine whether 
or not an interaction takes place. If an interaction does take place, a random number is 
generated and used to sample the inelasticity a, and the particle's momentum (pj before 
the interaction) is changed to p = Pi{l ~ a). Now the additional escape resulting from 
the momentum decreasing during the interaction is considered. The probability of the 
particle immediately escaping downstream due to "falling out of the box" is estimated. 



Prob. (escape, Pi — > p) = ^^'^^^ — 

L{Pi) 




-2, 



(20) 



and a random number is generated to determine if this happens. If the particle does 
escape, then w{t + At) particles with momenta p are binned in a histogram of accelerated 
particles, and a new particle is injected with momentum po and weight wq. If the particle 
does not escape, then the particle's momentum and weight are evolved through a new 
time step as described above. 

The resulting spectra are given in Figs. El and 0] give results for 5=1/3 and 1 for the 
cases for which we already have results for the corresponding continuous loss process, i.e. 
the same rioss(p) but with a — 0, for a = 0.005, 0.05 and 0.5, and (for simplicity) an 
inelasticity distribution which is constant for < a < 2a and zero elsewhere. In each 
figure I show spectra for combinations of r=1.75, 2 and 2.25, and rioss(p) oc p^^p^-^ and 
p^. The histograms show Monte Carlo results for the three inelasticity values, and the 
solid curve gives the analytic result for continuous losses which is found to be in excellent 
agreement with the d = 0.005 case. As is to be expected, where a strong pile-up (spike) 
is present for the case of continuous losses, it remains but is reduced if a = 0.05, and 
is completely removed if a = 0.5. In the last case, the spectrum generally continues 
to momenta above Pcut, and there can also be an enhancement (weak pile- up) at lower 
momenta. 

The spectra obtained above for the a = 0.5 and rioss(p) oc p° would be directly applica- 
ble to the case where proton acceleration is cut off as a result of hadronic interactions with 
matter above the threshold for multiple pion-production because the proton-nucleon in- 
elastic cross section increases only logarithmically with energy, and may be approximated 
as constant over, say, one decade in energy. Furthermore, for this case a not-unreasonable 
approximation for the inelasticity distribution is a uniform distribution between and 
1. Although other processes are normally assumed to cut off the spectrum, if the matter 
density is sufficiently high, it is in principle possible for hadronic collisions of protons 
with matter to play this role. Another possible application, again for the rioss(p) oc p^ 
case, would be for acceleration of electrons in relatively dense astrophysical matter in 
which relativistic bremsstrahlung could cut off the electron spectrum. Although a flat 
inelasticity distribution is not such a good approximation in this case, this process does 
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involve discontinuous energy losses comparable to the electron's energy, and so the results 
for d = 0.5 may serve as a guide. 

Perhaps, the most interesting application, again for the a = 0.5 and rioss(p) oc p° case 
would be for acceleration of high energy protons in a dense radiation field to energies 
well above the pion photoproduction threshold where rioss(p) is approximately constant 
and a ranges from 0.2 to 0.5 depending on proton energy. Near threshold rioss(p) is 
strongly dependent on momentum, and just above threshold the momentum dependence 
changes smoothly and passes through rioss(p) oc p and rioss(p) cx: p^/^ before becoming 
approximately constant well above threshold. For proton acceleration in a radiation field 
where pion photoproduction cuts off the spectrum, one can expect to see features of the 
type shown in Figs. OHU depending on the cut-off momentum, and momentum dependence 
of the diffusion coefficient. This process will be discussed in greater detail in the next two 
sections for the case of black-body radiation. 

5 Cut-off momentum for acceleration in black body 
radiation: sites of UHE CR 

Although I shall concentrate here on cut-offs due hadronic interactions with photons of 
the CMBR during acceleration of UHE CR protons, the results presented will be equally 
applicable to proton acceleration in any black body radiation field provided the proton 
momenta and interaction time-scales are scaled appropriately. For example, they may 
be applicable to proton acceleration in quasars in the nearly black-body accretion disk 
radiation. For an arbitrary temperature T = 2.725 T2. 7 K, the momenta quoted should be 
replaced by P/T2.7. 

In the case of acceleration of UHE CR, interactions with the CMBR may be important 
if the magnetic field is less than ~ 10^^ G. In this case, the gyroradii are rather large, 
and possible acceleration sites could be hot-spots in lobes of giant radio galaxies 12^] , 
or MHD generated electric fields along the whole length of an AGN jet such as that of 
Cen A 121] . Inevitably, acceleration times will be long, and one should also include redshift 
losses, and Bethe-Heitler pair production which has a somewhat lower threshold than 
pion photoproduction. In Bethe-Heitler pair production, the energy lost per interaction 
is small, less than 1%, and so this process can be treated adequately as a continuous loss 
process. The total momentum loss rate for these three processes is given in Fig. Efa) 
together with the momentum-gain rates required to reach various cut-off momenta for 
(5=1/3, 1/2 and 1. 

By plotting magnetic field vs. size of various astrophysical objects, Hillas [22] identified 
possible sites of acceleration of UHE CR based on whether or not the putative source could 
contain the gyroradius of the accelerated particles, and on the likely velocity of scattering 
centres in these sites (arguments similar to those that led to Eq. Q). Following Hillas |22] 
one finds that possible sites included neutron stars (10^-lO^^G), gamma ray bursts and 
active galactic nuclei (lO'^-lO^G), and lobes of giant radio galaxies and galaxy clusters 
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This identification of possible sources does not take account of energy losses 
(synchrotron) and interactions (Bethe-Heitler and pion photoproduction) which can cut 
off the spectrum, and so apply an additional constraint. In the same paper, Hillas |^ used 
this additional constraint to narrow the field of possible sources to radio galaxy lobes and 
galaxy clusters, and showed that to accelerate protons to ~ 10^°eV large regions (~Mpc) 
of relatively low magnetic field (less than ~ lO^^G) are needed, ruling out high magnetic 
field regions for the origin of UHE CR. For the higher magnetic fields in this range the 
spectrum is cut-off by synchrotron loss, and would be as discussed in Section 3. Thus, 
one of the very few plausible acceleration sites of UHE CR may be associated with the 
radio lobes of powerful radio galaxies, either in the hot spots [22i or possibly the cocoon 
or jet [SI- 

To estimate cut-off momenta (or energy), one needs plausible values for the accel- 
eration rate. The following values [2^ for the acceleration rate parameter in Eq. 2] are 
used: maximum possible acceleration rate ^(pcut)=l; plausible acceleration at perpendic- 
ular shock with speed 0.1c, ^(pcut) ~ 0.04, and plausible acceleration at parallel shock 
with speed 0.1c, ^{pcut) ~ 1.5x10"^. The latter value is consistent with that of Biermann 
and Strittmatter who considered DSA of protons in radiation and magnetic fields 
of active galactic nuclei. Assuming a Kolmogorov spectrum of turbulence, equation 6 of 
Biermann and Strittmatter [22] would give an acceleration rate parameter at the maxi- 
mum energy of ^ = 0.08/?^6 where (3sC is the shock velocity and b is the ratio of turbulent 
to ambient magnetic energy density; for /?s=0.1 this gives ^(pmax) ~ 8x10^^6. Based on 
the total momentum loss rate given in Fig. Efa), the proton cut-off momentum is plotted 
in Fig. El^b) as a function of magnetic field of protons for three adopted ^-values (chain 
lines are for constant Larmor radius as labelled). This plot, which is based on that in 
ref. j2S]; clearly shows that to accelerate protons to ~ lO^'^eV large regions of relatively 
low magnetic field ~ 10^''-10^^G are needed, apparently ruling out high magnetic field 
regions for the origin of UHE CR (see also ref. ^3)- ^^^^ that, in principle, protons 
can be accelerated up to ~ 5 x lO^^eV in Mpc scale region with ~ 10~^G. 

Fig. IS] is the "Hillas plot" with constraints added corresponding to the three curves 
in Fig. El^b), and where the chain lines give constant proton energy values as indicated; 
sources to the right of the solid curves being excluded. A possible exception to this is 
in the case of relativistically beamed sources (e.g. for AGN see ref. and for GRB 
see ref. jSO]) where neutrons emitted along the direction of relativistic motion can be 
Doppler boosted significantly in energy. Another possible exception is the case of so 
called "one-shot" mechanisms (e.g. [SIIE21) where a particle is accelerated by an electric 
field along a nearly straight path which is essentially parallel to the magnetic field such 
that curvature and synchrotron losses are negligible. Suggested sites include polarization 
electric fields arising in plasmoids injected into a neutron star's magnetosphere jHS] and 
magnetic reconnection in the magnetosphere of accretion induced collapse pulsars fM]. 
Another possibility is plasma wakefield acceleration, i.e. acceleration by collective plasma 
waves, possibly in the atmosphere of a GRB, or "surf-riding" in the approximately force- 
free fields of the relativistic wind of a newly born magnetar jSHl- In these cases it is 
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unclear whether the requirements of neghgible radiation losses can be met. I will next 
discuss in some detail shape of the spectrum, where the cut-off is due to photoproduction 
losses, for a large range of cut-off momenta. 



6 Spectrum for cut-ofF due to photoproduction loss 

Treating pion photoproduction as a continuous loss process, the escape rate and acceler- 
ation rate are given by 



'^esc(p) = (r - l)rioss(Pcut) 
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and the spectral index is 
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(23) 



(24) 



By examining the gradient of the solid curve in Fig. Efa) one sees that q{p) has a strong 
peak at ~ 10^° eV, and this will give rise to a strong pile-up in the spectrum if the 
cut-off momentum is close to ~ 10^'' eV/c. The resulting spectra for the continuous 
loss approximation are shown in Fig. d and do indeed show this strong pile-up if pcut ~ 
10^° eV/c. Pile-ups also occur near other pcut values, particularly if pmax ^ Pcut, and are 
strongest if 5 = 1. However, when the constant escape term (corresponding to, e.g., escape 
due to finite size of accelerator) dominates escape near pcut, then the resulting steepening 
in the spectrum is the prominent feature in all cases (except where pcut ~ 10^° eV/c when 
there is always a strong pile-up). Of course, the above treatment assumes Bethe-Heitler 
pair production and pion photoproduction to be continuous loss processes. While this is 
a reasonable approximation for Bethe-Heitler pair production, it certainly is not for pion 
photoproduction, and in this case the sharp cut-off and pile-up will be partly smeared. 
Monte Carlo simulations have been performed to see to what extent any pile-ups shown 
in Fig. m are washed out. 

In the Monte Carlo simulations, redshifting and Bethe-Heitler pair production were 
treated as continuous loss processes, but hadronic collisions of protons and neutrons were 
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treated by the Monte Carlo method using the SOPHIA event generator jSH]. Neutrons 
produced, e.g. in jrj nvr"*", may decay inside the box into protons which continue 
to be accelerated, or they may escape from the box to decay outside into cosmic rays 
depending on the neutron's Lorentz factor and the physical dimensions of the acceleration 
region. Since the simulations are for a given pcut, Pmax and 5, and the size of the box 
depends on the (unknown) normalization of the diffusion coefficients, I have performed 
simulations separately for the two extreme cases. In the first case the cosmic ray injection 
spectrum is made up only of protons escaping from the box, whereas in the second case 
the cosmic ray injection spectrum is made up of all neutrons produced inside the box 
plus protons escaping from the box. High energy gamma rays and neutrinos also result 
from photoproduction during acceleration (see ref. |SZ1 for acceleration of protons and 
production of neutrinos during acceleration in ~10^ K blackbody radiation). 

The Monte Carlo results for proton acceleration in a black body radiation field with 
r = 2 and T2,iPmaxC = lO^^eV are shown in Fig.|Sl with the left column being for neutrons 
decaying inside the box, and the right column being for neutrons escaping from the box 
before decaying. Upper thin curves are for escaping protons, and lower thin curves are for 
neutrinos (all flavours) produced during acceleration. Comparing the T2,7Pms.xC = lO^^eV 
case (right column in Fig. [7j) with the case of neutron decay inside the box (left column in 
Fig. IHI) confirms the presence of the general features expected from treating interactions 
as a continuous loss process. The sharp features of Fig. [3 are somewhat rounded, and 
the spectrum contains particles with momenta exceeding pcut, but this is to be expected 
because of the distributions of inelasticity and interaction time. 

If neutrons escape from the box before decaying, the protons from neutron decay will 
not be further accelerated and this explains why the spectra in this case are lower than the 
case of neutrons decaying inside the box. In some models, protons are trapped outside of 
the acceleration region but inside the astrophysical source, and so it is of interest to know 
separately the spectrum of escaping neutrons produced during the acceleration process 
(there will also be be a component due to pion photoproduction of the trapped accelerated 
protons). The contribution of escaping neutrons to the total escaping proton flux is show 
in Fig. ISl by the thick curves. 

For the first two spectra, with PcutC = 1-3 x lO^^eV, the cut-off is due to Bethe- 
Heitler pair production, and so there is no significant neutrino or neutron component. If 
PcutC = 1-3 X lO^^eV, there will be a strong pile-up. For PcutC > 3 x 10^°eV, the nature 
of the spectrum above 3 x lO^'^eV depends strongly on the momentum dependence of 
the diffusion coefficient. For k oc p^^^, and neutrons escaping, the spectrum essentially 
breaks from ~ below the pile-up feature, to ~ p~^-p~^ above the pile-up feature 
(corresponding to PcutC = 10^^-lO^^eV). If neutrons decay inside the box, the flux above 
the pile-up feature is higher and more curved. Nevertheless, for PcutC > 3 x 10^°eV and 
with K oc p^^^ the flux at p ~ pcut is substantially lower than expected without interactions. 
In other words, the effect of interactions cuts in at much lower momenta than p^nt, and 
this is most pronounced for k oc p^^^. For k (x p the spectra above 3 x lO^^eV generally 
continue as roughly the original ~ power law and suffer an almost exponential cut-off 
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at p ~ Pcutj with the cut-off occurring at a shghtly lower momentum in the case of neutron 
escape. As would be expected, the spectra for k oc p^/^ are intermediate between the the 
p^/^ and cases. 

6.1 Spectrum after propagation 

Unless the source is very near by, and the propagation time is less than ~10 Mpc/c, then 
the observed spectrum will be severely affected by propagation. This is illustrated in Fig.lHl 
which shows spectra of protons and neutrinos (all flavours) escaping from the acceleration 
region (dotted curves) and after propagation for time 100 Mpc/c for the cases given in 
Fig. ISl Results (not shown) for 50 Mpc/c are only marginally different, implying that 
cascading is essentially complete by 100 Mpc/c. The proton spectra after propagation 
for all cases are hardly distinguishable, as they all show, as expected, a prominent GZK 
cut-off at 10^°eV with a pile-up just below the acceleration cut-off. Hence, it appears that 
one can conclude little about the acceleration cut-off energy and hence physical conditions 
of the acceleration region, by directly observing the spectrum of UHE CR. 

Information is preserved in the neutrino spectrum, to some extent, about the proton 
cut-off energy and the momentum dependence of the diffusion coefficient during accel- 
eration. During propagation, the energy of protons above the GZK cut-off is converted 
into secondary particles, including neutrinos which propagate essentially unhindered. As 
can be seen, the cut-off in the neutrino spectrum after propagation is directly related to 
the cut-off in the accelerated proton spectrum, and the spectral shape retains some infor- 
mation about the shape of the cut-off of the accelerated proton spectrum, and hence the 
physical conditions of the accelerator. The spectra plotted are effectively E'^dN/dE and 
so are proportional to the energy flux of particles per logarithmic energy interval, or the 
"spectral energy distribution" (SED). In Fig. irUTa) I plot the energy at the peak of neu- 
trino SED after propagation vs. proton acceleration cut-off momentum, and in Fig. ITUT b) 
the maximum neutrino energy, which I define here as the energy at which the neutrino 
SED is 1/e of its peak value, vs. proton acceleration cut-off momentum. As can be seen, 
the maximum neutrino energy is a more useful measurement for determining the proton 
acceleration cut-off momentum as it is less subject to statistical fluctuations (in real or 
Monte Calo data) and is more sensitive to pcut- 

The inferred proton acceleration cut-off momentum is sensitive to whether or not the 
neutrons escape from the acceleration region, particularly for the case of a Kolmogorov 
spectrum of turbulence. However, one can still use Fig. ITUT b) to put limits on the proton 
cut-off energy. For example, the observation of UHE neutrinos with E^ ^^^y^ = 3 x lO^'^eV 
would imply a proton cut-off energy of 3-4xl0^^eV (5=1), 5x 10^^-lO^^eV (5=1/2), and 
1.5xl022-1023eV (5=1/3). 
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7 Conclusion 



In this paper I have discussed the shape of the high-energy part of the spectrum of 
particles accelerated by DSA in the presence of (i) additional diffusive escape from the 
accelerator, (ii) continuous energy losses, (iii) energy changes arising from interactions. 
I find that the shapes of spectral features (roll-offs, pile-ups, cut-offs) depend strongly 
on the momentum-dependence of the acceleration rate, which in turn depends on the 
momentum-dependence of the diffusion coefficient, k oc ^ and on the physical process 
responsible for inhibiting acceleration. Hence, the shape of the high-energy part of the 
spectrum is important because it can give clues to physical conditions at the acceleration 
site. In some situations the accelerated spectrum could, in principle, be observed directly, 
but more usually the particle spectrum (if observed) would be affected by propagation 
from source to observer. Nevertheless, some features of the spectrum after acceleration 
could remain, or the spectrum of radiation emitted by the accelerated particles would 
retain some sensitivity to the features of the accelerated particle spectrum. 

In the case of continuous losses where the spectrum is cut off at the momentum at 
which acceleration rate equals loss rate, pile-ups occur for p^-losses for F < 2, and for 
r=2.5 if 5 = 1. However, for p^-losses, there is a fairly strong smooth cut-off except 
for the flattest spectrum considered, r=1.5 with 5 = 1. In the case of additional escape 
processes inhibiting acceleration at momenta comparable to or below the cut-off momen- 
tum, no pile-ups are seen, and in all cases the spectra roll off smoothly at the momentum 
where the additional escape becomes comparable with the normal escape downstream. 
One application of this could be in detailed fitting the X-ray spectra of SNR ^38j, or 
other sources in which synchrotron radiation is observed. Another example where these 
subtleties could be important is in acceleration in AGN jets where, for some physical con- 
ditions, the spectrum could be cut-off either by adiabatic losses as the jet expands, or by 
leakage out of the jet. In the jet-frame a reasonable order of magnitude approximation for 
the time scales of both processes is the jet radius divided by c. Depending on which pro- 
cess dominates, and the physical parameters, the spectrum could be strongly cut off with 
spectrum affected well below the nominal cut-off momentum (adiabatic losses), or have 
a more gentle decline staring just before the nominal cut-off momentum and extending a 
bit above it (diffusive escape). 

In the case of interactions, the fate of the spectral features depends largely on the 
mean inelasticity of the interaction. For very small inelasticity, the results are close to 
the equivalent continuous loss case, but if a ~ 1 pile-ups and cut-offs are smoothed and 
particle spectra extend beyond the nominal cut-off. I have considered in detail the case 
of pion photoproduction by protons on black body radiation, focusing on acceleration of 
UHE CR in the presence of the CMBR. Using plausible acceleration rate parameters, I 
have re-examined the physical conditions necessary to accelerate protons to UHE and, 
by considering a "Hillas plot" with constraints from photoproduction and synchrotron 
radiation added, confirm that the most likely of the originally proposed acceleration sites 
are in hot-spots in lobes of giant radio galaxies, as inferred by Hillas 22 and others 
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subsequently. However, relativistically moving sources such as AGN jets and GRB do not 
need to accelerate particles directly to UHE as Doppler boosting from the plasma frame, 
in which the acceleration occurs, to the host-galaxy frame can increase the energies of 
escaping particles, for example, by ~30 for AGN jets and ~1000 for GRB. Thus, even 
though these sources appear to be ruled out on the modified Hillas plot, they could well 
be viable. 

If the UHE CR are indeed accelerated, and do not result from a top-down scenario, 
then the high energy end of the accelerated particle spectrum could contain clues of the 
physical conditions of the acceleration region. The spectra on acceleration for different 
physical conditions are quite distinct. If the source of the UHE CR is Galactic, e.g. a 
GRB which occurred in our Galaxy in recent pre-historic times (say within a few tens of 
thousands of years), then one might directly observe pile- up features, possibly at super- 
GZK energies. If the sources are extragalactic the particle spectrum one would observe 
would almost certainly have had all information about the part of the spectrum above 
~ 10^°eV, i.e. beyond the GZK cut-off, removed. However, UHE neutrinos produced 
during propagation can be used to estimate at least the cut-off energy of the protons 
just after acceleration, if some reasonable assumptions about momentum dependence of 
diffusion coefficients and size of the accelerator can be made. The neutrinos produced 
during propagation will not point back directly to the acceleration region, but will appear 
to come from a halo of size somewhere between the proton gyroradius in intergalactic 
space around the source and the pion photoproduction loss distance,i.e. ~1 to 15 Mpc, 
depending on the magnetic field environment around the source. At the 16 Mpc distance 
of M87 this would mean neutrinos would arrive from this source at angles up to 4°- 
90°, whereas for Cen A at only ~4 Mpc, one might directly observe cut-off features in 
the UHE CR spectrum, but one would observe only those neutrinos from the direction 
of Cen A produced during acceleration unless the accelerated cosmic rays are efficiently 
trapped near to Cen A. For more distant source candidates, it might be possible to 
identify some currently active sources of UHE CR. Finally, I mention that if heavy nuclei 
are accelerated then photodisintegration may take place during acceleration and as well 
as during subsequent propagation. Simulations are planned for this case and will be the 
subject of a future paper. 
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log(x) log(a;) 

Figure 1: Differential energy spectra for xq <^ 1 <^ x^ax ior £2 — 1/2, 6 as indicated, and 
r = 1.5 (upper curves), 2.0 (middle curves) and 2.5 (lower curves). Results are shown for 
(a) dp/dt oc p^, (b) dp/dt oc p. 
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Figure 2: Differential energy energy spectrum for the case of F = 2 for (a) S = 1/3, (b) 
5 — 1/2, (c) S — 1, and for (d) F = 2.2 and 5 — 1/3. In each case £2 — 1/2, and results are 
shown for momentum-dependent escape at Pmax and Pcut/Pmax = 10~^ (leftmost curve), 
10""^, . . . , 10^ (rightmost curve). 
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Figure 3: Spectra for 5 = 1/3, £2 = 0.5 and combinations of r=1.75, 2 and 2.25, and 
^'ioss(p) oc and as indicated. In %9,ch case the histograms show Monte Carlo 

results for q;=0.005 (thin histogram), q;=0.05 (intermediate histogram) and a— 0.5 (thick 
histogram), and the solid curve gives the analytic result for continuous losses. 
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Figure 4: As Fi^;. El except 5 = 1. 
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Figure 5: (a) Energy loss rate of protons rioss(p) in black body radiation of temperature 
T = 2.725 T2, 7 K with losses by Bethe-Heitler pair production, pion photoproduction and 
redshifting ^21- Note that redshift losses are included to illustrate their importance in the 
case of UHE CR acceleration in the 2.725 K CMBR. For other temperatures, the rate for 
Bethe-Heitler pair production plus pion photoproduction (chain curve) should be used; 
the dot- dot- dot- dash curve gives the loss rate for pion photoproduction only. Acceleration 
rates needed to reach cut-off momenta Pcutc/T2.7 = 10^^, 10^^-^, . . . , 10^^ eV are shown by 
the thin lines: short dashed {6 = 1/3), long dashed {S = 1/2) and solid {6 = 1). Dotted 
lines show synchrotron loss rates for B/T2J = 3 x lO^^G and 10~^G. (b) Maximum energy 
as a function of magnetic field of protons for maximum possible acceleration rate ^ = 1 
(upper solid curve), plausible acceleration at perpendicular shock ^ = 0.04 (middle solid 
curve), and plausible acceleration at parallel shock 1^ = 1.5 x 10~^ (lower solid curve). 
Dashed curves are limits from Bethe-Heitler pair production and pion photoproduction 
only (solid curves include synchrotron loss). Dot-dot-dot-dash curves are lines of constant 
Larmor radius as labelled. 
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Figure 6: "Hillas plot" showing (chain curves) magnetic field vs. gyroradius for proton 
momenta 10^^, 10^^, . . . , lO^^eV/c. The solid curves correspond to those in part Fig.E^b) 
such that the parameter space of accelerated particles is to the left of the curve corre- 
sponding to the chosen acceleration rate parameter. Typical size and magnetic field of 
possible acceleration sites (taken from Hillas are shown for neutron stars (ns), white 
dwarfs (wd), sunspots (ss), magnetic stars (ms), AGN (ag), interstellar space (is), su- 
pernova remnants (sn), radio galaxy lobes (rg), galactic disk (d) and halo (h), clusters 
of galaxies (cl) and intergalactic medium (ig). Typical jet-frame parameters of the syn- 
chrotron proton blazar model [21] and gamma ray burst model jHOl are indicated by the 
open squares labelled "bl" and and "gb". 
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Figure 7: Differential momentum spectra of protons accelerated in the CMBR with con- 
tinuous losses by redshifting, Bethe-Heitler pair production and pion photoproduction for 
r = 2, and for 6 and PmaxC as indicated. In each case results are shown for Pcutc = 10^^ 
(leftmost curve), 10^^-^, . . . , 10^^ eV (rightmost curve). 
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6=0.33, neutrons escape box 




6=1.00, neutrons decay in box 6=1.00, neutrons escape box 
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Figure 8: Differential momentum spectra of protons accelerated in the CMBR with con- 
tinuous losses by Bethe-Heitler pair production and redshifting, and pion photoproduction 
interactions simulated using the SOPHIA event generator for F = 2, PmaxC = 10^^ eV, 
5 as indicated, and for neutron decay inside the box (left column) and neutron escape 
from the box (right column). Upper thin curves show the proton flux for PcutC = 10^^ 
(leftmost curve), 10^^'^, . . . , 10^^ eV (rightmost curve). Lower thin curves total neutrino 
flux (all flavours), and thick curves show spectra of escaping neutrons, in both cases for 
PcutC = 10^° (leftmost curve), 10^°-^, . . . , 10^^ eV (rightmost curve). 
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I (5 = 0.33, neutrons decay inside box:: 6=0.33, neutrons escape from box: 
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6 = 0.50, neutrons decay inside box: 
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6=0.50, neutrons escape from box: 




: 6=1.00, neutrons decay inside box:: 6=1.00, neutrons escape from box: 
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Figure 9: Spectra of protons and neutrinos (all flavours) escaping from the acceleration 
region (dotted curves) and after propagation for time 100 Mpc/c (solid curves) for PcutC = 
10^° (leftmost curves), 10^°'^, . . . , 10^^ eV (rightmost curves) for the cases given in Fig. |H| 
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Figure 10: (a) Energy at peak of neutrino SED after propagation for time 100 Mpc/c vs. 
proton cut-off momentum, and (b) maximum neutrino energy defined as the energy at 
which the neutrino SED, after propagation for time 100 Mpc/c, is 1/e of its peak value vs. 
proton cut-off momentum. Thin curves are for neutron decay inside acceleration region, 
and thick curves are for neutron escape from acceleration region; results are shown for 
5=1/3 (solid curves), 1/2 (dotted curves) and 1 (dashed curves). 
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